Thirring solitons in the presence of dispersion 



Alan R. Champneys 
Department of Engineering Mathematics, 
University of Bristol, Bristol BS8 1 TR UK 



Boris A. Malomed 
Department of Interdisciplinary Studies, 
Faculty of Engineering, Tel Aviv University, Tel Aviv 69978, Israel 



Mark J. Friedman 
Department of Mathematics, University of Alabama,IIuntsville, Alabama 
(DRAFT February 2, 2008) 



Abstract 



The effect of dispersion or diffraction on zero-velocity solitons is studied for 
the generalized massive Thirring model describing a nonlinear optical fiber 
with grating or parallel-coupled planar waveguides with misaligned axes. The 
Thirring solitons existing at zero dispersion/diffraction are shown numerically 
to be separated by a finite gap from three isolated soliton branches. Inside the 
gap, there is an infinity of multi-soliton branches. Thus, the Thirring solitons 
are structurally unstable. In another parameter region (far from the Thirring 
limit), solitons exist everywhere. 
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The massive Thirring model (MTM) |T|] is a completely integrable Lorentz-invariant 
model of classical field theory, which supports exact soliton and multi-soliton solutions 
A generalization of the MTM, which we call the optical model (OM), was introduced in 
to describe interaction between right- and left-traveling waves in a nonlinear optical 
fiber with a grating. Resonant Bragg scattering and cross-phase modulation (CPM) give 
rise, respectively, to linear and nonlinear couplings between the two waves. OM additionally 
includes self-phase-modulation (SPM) nonlinear terms, making it nonintegrable and destroy- 
ing Lorentz invariance. Nevertheless, a family of exact one-solitons can be found with 
arbitrary velocity and internal amplitude ( "soliton" hereafter means solitary wave, and "n- 
soliton" is one with n peaks). Recently, Bragg solitons have been observed experimentally 
in a fiber with grating [^]. 

Both MTM and OM neglect dispersion of the medium, solitons being supported by an 
effective dispersion induced by the linear coupling. In physical media, however, material 
dispersion is present. The aim of this work is to examine the influence of such dispersion D 
on the Thirring solitons (TS). This first study treats only zero- velocity solitons, which are 
essentially the same in MTM and OM. Results for finite- velocity ( walking P) sohtons, to 
be presented elsewhere, are more complicated technically but not drastically different (see 
below). The zero- velocity solitons are most intriguing physically, as they imply complete 
dynamical self-trapping of light on the grating. We will conclude that TS are structurally 
unstable (which does not imply their dynamical instability), being separated by a finite gap 
from the nearest branch of fundamental-solitons for D > 0, and with no solitary waves at 
all for D < 0. Within the gap, we find infinite sequences of two-solitons that are bound 
states (BS's) of the fundamental ones. Although likely to be dynamically unstable |^, BS's 
are worth studying to delimit the existence domain of fundamental solitons, see Figs. 2,4 
below. 

The generalized MTM including dispersion terms is 

iut + iux + Duxx + (c^l^^P + I'i^r) M + f = 0, (1) 

ivt — ivx + Dvxx + {\u\^ + f + M = 0, (2) 

where u{x,t) and v{x,t) are the complex amplitudes of the counterpropagating waves, x 
and t are the coordinate and time, D is the coefficient of spatial dispersion, and a is the 
relative SPM coefficient, which is zero for MTM, and | for the OM case. Besides fibers 
with grating, the model (|I]), (H) can be applied to stationary fields in two parallel tunnel- 
coupled planar nonlinear waveguides. In that case t and x are the propagation distance and 
the transverse coordinate respectively, the terms ±iux account for misalignment of optical 
axes in the two cores, D is an effective diffraction (not dispersion) coefficient, and the CPM 
terms must be omitted (see, e.g., 0). Actually, the latter realization of the model is closer to 
experiment, as optical axes misalignment is a powerful control parameter enabling rescaling 
of physically realistic systems into the form (0), (0) pj. In contrast, for fibers with grating, 
a simple estimate shows that dispersion may not be conspicuous unless the spatial width of 
the soliton is comparable to the grating period, i.e. the wavelength of light, when Eqs. (1) 
and (2) are not applicable [^] 

Essentially the same model governs interaction of two circular polarizations of light in a 
nonlinear fiber, in which the linear coupling is induced by the birefringence, and the group- 
velocity difference by a fiber's twist (see the review [|10|)- In untwisted fibers, interaction 
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between linear polarizations is described by similar models but with linear coupling replaced 
by a cubic four- wave-mixing term (see |11] for a family of walking solitons in the latter 
model). Our approach is different; instead of starting from solitons of decoupled nonlinear- 
Schrodinger equations with the couplings treated as perturbations, we start from the TS of 
the strongly coupled system with dispersion or diffraction being a perturbation. 

Being interested here only in the zero-velocity solitons, we substitute into Eqs. (|I]) and 
(g) u{x,t) = e-*'^*f/(x), V = e-''^^V{x) to obtain the coupled ODEs 

DU" + lU' + ujU+ ((T|f/p + \V\'^)U + V = (3) 
DV" - iV + UJV + {\U\^ + a\V\^)V + U = 0, (4) 

the prime standing for d/dx. In this notation, the TS occur at D = and \uj\ < 1. Eqs. 
and (^) are equivalent to an 8th-order dynamical system with two integrals of motion: the 
Hamiltonian 

H = D(|?7f + ) + + + (a/2)(|[/|^ + |V^|^) 

+ \U\^\V\^ + {UV* + VU*), (5) 

and the "angular momentum" , generated by invariance with respect to the continuous phase 
transformation, 

M = D (UU'* - U*U' + W* - V*V') + \V\^ - \U\^. (6) 

This Hamiltonian system has several discrete symmetries: the odd symmetry Z : (f/, V) — > 
(— [/, —V), two other Z2 ones Zi : f/ ^ V* , Z2 :U —V*, and four reversibilities 

R ■ {U, U', V, V) ^ {U\ -U'\ V\ -V'*), -x^-x, (7) 
S : [U, U') ^ [V, -V), : X -> -X , (8) 

along with their odd images ZR and ZS. 

The first step in locating solitary waves is to solve the linearized problem, assuming 
solutions ~ e'^^ . This problem, solved exactly, gives a set of double eigenvalues: 

(d^a^ + 2Du\^ + + ^2 _ ^ . (9) 

Eq. (^) defines four regions on the plane {D,uj} with different types of eigenvalues (see Fig. 
1). Solitary- waves with exponentially decaying tails are only possible in regions I, II, and 
III (and their images for D < 0), where eigenvalues with nonzero real part occur. 

We notice that Eqs. (^ and (^) are compatible with the reduction U = V* . This results 
in a single equation for U{x), 

DU" + lU' + ujU+{1 + (y)\U\^U + f/* = 0, (10) 

equivalent to a real fourth-order ODE system. All the zero-velocity solitons in MTM and 
OM obey exactly the same reduction, and a simple argument based on consideration of the 
unstable manifolds shows that all possible zero- velocity solitons to ([I|), (0) within region II 
are trivially related to solutions of ([T0|) by rotation in the (f/, V) plane. Henceforth, we set 
(7 = because a can be scaled out from Eq. (PD]). Furthermore, for Eq. (|T^), S = R, and 
the "angular momentum" (^ identically vanishes. The eigenvalues of the corresponding 
linearized equation are given by Eq. (P), but are all single, i.e.. Fig. 1 remains fully relevant. 
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FIG. 1. Parameter regions for D > with different types of eigenvalues of the hnearized Eqs. 
@ and as ihustrated by the insets. The curve delimiting region II is -D + 1/4D — uj = 0. At 
the point D = ^, the curve is tangent to the horizontal uj = —1. The picture for D < is obtained 
by rotating the figure by 180°. 

The soliton is a homoclinic-to-zero solution to Eq. (plQl). According to general theorems 
1^ , in region III, where [/ = is a saddle-center fixed point, homoclinic trajectories that 



are symmetric under a reversibility are of codimension one (non-symmetric homoclinic tra- 
jectories are of codimension two). Hence solutions can exist only on isolated curves in the 
{D,uj} parameter plane, the number of which may be finite or infinite. Moreover, given a 
sign condition on the quadratic part of the Hamiltonian, each curve will be accompanied 
by an infinite accumulation of curves on which BS's exist |]T^. In contrast, in regions I 
and II, where the fixed point f/ = is hyperbolic, homoclinic trajectories are generic, i.e., 
they occur uniformly in two-dimensional parameter regions [0. But region III is of most 
interest, as it abuts the segment {D = 0, |ti;| < 1} on which the TS solitons exist. 

To obtain solutions, we use robust numerical methods for solving two-point boundary- 
value problems on a truncation of an infinite x-interval with boundary conditions placing the 
solution in the stable or unstable eigenspaces at the origin; see |]14| and references therein. 
Continuation of solutions with respect to parameters is carried out using the software auto 
|T5|, specifically exploit the reversible structure of (0). 

Our main findings are summarized in Fig. 2. Here, three solid curves represent the 
isolated loci of fundamental or primary (single-humped, in one component) solitons, and 
the dashed curves are a small sample of loci of their two-humped BS's. All primary-solitons 
are reversible with respect to the transformation ZR, see Eq. (0) ; we have found no evidence 
of any i?-reversible solutions. In panel (b), we use, instead of the frequency a;, the soliton's 
energy E = \ U{x)\'^dx. Typical examples of one-solitons are displayed in Fig. 3. and 
typical two-humped BS's are shown in Fig. 4 (only half of each two-soliton is shown in this 
figure). 
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FIG. 2. A two-parameter bifurcation diagram for fundamental solitons (solid curves) and 
two-soliton bound states (dashed curves) on the planes {D,uj) (a) and (D,energy) (b). 
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FIG. 3. The fundamental solitons at the points of intersection of the primary-soliton branches 
(Fig. 2) by the line uj = -0.8. 

Each of the primary branches in Fig. 2 (labeled 1-3) appears to bifurcate at zero soliton 
amplitude from the line u = 1, although there are numerical difficulties in computing right 
up to this singular limit. The D -values of these three bifurcations at uj = 1 are D = 0.50, 
0.20 and 0.11 to two decimal places. A straightforward calculation of the sign condition 
in [T^ on the Hamiltonian (^ implies that curves of Zi?-reversible BS's must accumulate 
on each of primary curves from both sides (e.g. we have found BS branches 9-12 and 4,5 
accumulating on branch 1 from the right and left respectively). Branches 9-12 are also part 
of a larger sequence we have computed which for fixed D accumulates on u = 1. 

Three-solitons and higher-order BS of the primary solitons can also be found, the three- 
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humped ones accumulating on two-solitons, etc. in accord with the theory [jT2[. We do not 
describe these objects because it is unhkely that even the two-sohtons may be dynamically 
stable in the PDE system (|1|),(|^), while stability of the primary solitons is quite feasible 
0. However, stability analysis is deferred to another work. Homoclinic solutions were also 
sought for D < and |a;| < 1, but no evidence of primary or multi-humped ones was found. 

Looking at Fig. 2, there remains the crucial question whether there are any more primary 
branches to the left of that labeled 3. A seemingly plausible conjecture is that there is a 
self-similar structure of primary branches as one moves to the left in Fig. 2, i.e., infinitely 
many branches accumulating on the TS segment {D = 0, Ic^l < 1}, the branches 1,2 and 3 
being but the first three in the structure (note that, at least for D < 2, there cannot be any 
further primary solutions to the right of branch 1, because here the two-soliton curves 9-12 
form a barrier for them). However, careful numerical scanning of the parametric plane of 
Fig. 2 to the left of branch 3 has strongly indicated that the above hypothesis is false, in 
region HI there being no primary branches other than 1, 2, and 3. For u sufficiently close 
to —1, this assertion is substantiated as follows 

Fig 4. reports the results of a thorough numerical investigation of other possible solution 
branches at = —0.99, varying D between and |. We find that, to the left of branch 
3, an infinite sequence of mu/iz-soliton BS's occurs. Even though, because of numerical 
problems in the singular limit, we have only computed the corresponding two-solitons down 
to ~ 0.2, Fig. 4 clearly suggests accumulation of the sequence as D —>■ 0. The energy of 
the two-solitons remains finite, while the separation between the two bound pulses diverges 
~l/DasD^O (which explains the existence of TS in the limit D = 0). Thus, what does 
accumulate on the TS manifold at D — is an infinite sequence of multi-soliton branches, 
with no fundamental-soliton branch closer to the TS manifold than the branch 3 in Fig. 2. 

To support this numerical finding with qualitative arguments, consider what happens 
to the primary branches as they cross the line uj = —1 from above. For D > | this is a 
"harmless" transition, because the real eigenvalues of the linearized equations, that govern 
the decay of the homoclinic solution at |a;| —>■ oo, behave smoothly and they are bounded 
away from zero. A well-defined primary branch safely crosses u = —1 in this case, which 
for u < —1 describes a curve of "orbit-fiip" bifurcations (cf. [0). However, for D < ^, the 
corresponding eigenvalues vanish as ^ — 1, hence no smooth transition can take place. 
Thus, there may be no primary-soliton branches atO<l-|-u;<^l, D < ^. 
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FIG. 4. A sequence of two-soliton solutions at D 
shape of the solutions. 
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The most important result of this work is that there is a finite gap separating TS, exist- 
ing in the singular limit D = 0, from new solitons at -D 7^ 0. Thus, the Thirring solitons 
are structurally unstable against adding the dispersion or diffraction. A natural question 
is if there is a gap for solitons at a finite velocity. Preliminary numerical results give a 
positive answer, which is further supported by an argument that solutions to the ODEs 
describing the soliton's shape continuously depend on the parameters (including velocity), 
except at a singular point. The addition of dispersion to MTM is, obviously, a singular per- 
turbation because it doubles the system's order; however, nonzero velocity is not a singular 
perturbation. 

We mention finally results for solitons in regions I and II. As stated, here homoclinic 
solutions are generic, and a primary soliton branch can be path-followed continuously for 
all uj and D inside regions I and II. Inside region II it develops oscillations in its tails due 
to the complex eigenvalues. At the boundaries between regions I and III and II and IV, the 
solution disappears through a zero-amplitude bifurcation, as predicted by the appropriate 
normal-form analysis fl^. Other primary-soliton solutions have more complicated bifurca- 
tion diagrams; in both regions I and II, two- and multi-soliton BS's also occur. A detailed 
description of the complete bifurcation structure will be given elsewhere. 



Since the original submission of this paper, we have become aware of the preprint JIB 
containing new results on the dynamical stahiliiy of the solitons in OM without the dispersion 
terms. They demonstrate that, except for the integrable Thirring model case, all the solitons 
are subject to an instability which is too weak to have been observed in earlier numerical 
simulations. Note that a similar instability mechanism for solitons of OM was predicted 
non-rigourously in using a variational approximation. A dynamical stability analysis 
for the new solitons in the presence of dispersion found in the present work will be presented 
elsewhere. 

We appreciate valuable discussions with Y.S. Kivshar, G.G. Luther and D.E. Pelinovsky. 
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